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I. Introduction 

The design of adaptive control systems based upon the application of 
Liapunov theory has mainly been concerned with the idealized situtation in 
which the system is free of disturbance [1-3]. The effect of parameter devia- 
tions for the case in which the number of adaptive gains is less than the full 
set required for complete adaption has also been ignored. Based on these 
simplifications it can be shown that the tracking error of the adaptive system 
is asymptotically stable. 

In recent years efforts have been made to modify the Liapunov design in the 
interests of practicality and generalization. The main results have been to 
reduce the number of derivatives of the output variable which need be measured, 
and to generalize the synthesis procedure to permit m.ore rapid convergence of 
the tracking error [4]. At least one effort has been made to apply the 
design to a practical problem [5]. However before the design can be considered 
to have real engineering significance, the effects of incomplete zciaptibh 
disturbance and saturation must duly be considered. 

The purpose of this paper is to consider the first two of these items 
with respect to the particular adaptive configuration shown in Figure 1, In 
order to obtain results which can be readily Interpreted, the plant was assumed 
to be of less general form than is required by the existing theory. However 
it is evident from the results obtained that stability problems nay manifest 
themselves in the presence of disturbance and adaption errors, and that these 
problems should not be ignored. 
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First a sufficient condition for boundedness of the tracking error in 
the presence of disturbance and adaption errors will be derived. From this 
result it is shown that boundedness of the tracking error nay not in itself 
quarantee boundedness of the adaptive gain parameters. Hence an independent 
analysis is required to ascertain the effect of disturbance upon stability 
of the adaptive gains. It is shown that the disturbance and the input signal 
can be so related that the adaptive parameters will in fact be unbounded. 

A simulation study is carried out with respect to a third-order plant. 

It is shown in a practical situation that incomplete adaption can lead to a 
reasonable result, as predicted by the derived error bound. It is further 
demonstrated that Instability due to the action of a disturbance can be brought 
under control if the input to the system is properly chosen, and that the 
frequency of the input signal has a significant effect upon the tracking errors. 

II* Description of the Adaptive Control System 

The tine-invariant linear plant to be considered is defined according to 
the state equation 

x = Ax + bu + c, r + d (1) 

- pr- -p -p - 

wherein the state vector x = [x. ] is of dimension n. A = [a,,], b = [b. ], 

— i p ij -p i 

^ = [c^] contain constant unknown parameters, and id = [d^^] is an unknown 
bounded disturbance. In assuming that the plant has but one (scalar) input, 
the possibility of adjusting the plant coefficients directly is 

ruled out. Hence, as will be seen, adaption is to be obtained by the use of 
adjustable gains whose outputs act through the control input u so as to cause 
the plant output to track that of a model. As a consequence of this asstimption 
it can be shown [6] that the state variables in (1) must be chosen as phase 
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var tables, i.e., x = x-, x, = .... k . = x , in which case the following 

1 z / j n— 1 • n ^ 

forms for A , b and d are obtained: 



In the ensuing analysis only the sign of must be known. 

The stable tine-invari^t model, which is assumed to have the same struc- 

( 2 ' 

ture as the plant, is defined by, ' • 

i " V" V 

where y = [y . ] . A = b = [3. ] have a form simular to A , b above, 

^ m ij ’ -m i P “P 

and r is a scalar input signal. If the error vector is now defined by 


e = X - X. 

then the error differential equation can be cast in the form 

e = A e + f 
— n— — 

where 


f=(A -A)x+(b -c)r-bu-d=Ax+6r-bu-d. 
— n P “ — p — p — — — — 


(3) 

(4) 


Because of the phase- variable assumption, it can be seen that A, 6 have the 



Therefore ^ has only one non-zero element, namely 
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With reference to Figure 1 the objective is now to design an adaptive 
controller so as to realize an asyiaptctic bound on the nom of the tracking 
error, 1 1 e| | , where £ is goverened by (3). Ideally it is desired that 11^11 
should go to zero. In general it is seen that this may require n+1 adaptive 
gains, however, as will be seen in a case example, such an extravagance may 
not always be justified because deviation in certain of the plant parameters 
can have relatively small effect upon the system response, and each adaptive 
gain increases the level of measurement noise entering the system. The first 
objective is then to determine a bound on the tracking error when less than 
n+1 adaptive gains are incorporated in the system, and then to show that a 
disturbance id entering the system can cause instability. 

Ill, Description of the Control Law 

According to Parks [1], the design objective will be realized by 
synthesizing a Liapunov function. Thus starting with the scalar function 

V = e''^ P e + (5) 

in which P = to be positive definite symmetric, and ^ is a 

parameter dependent vector to be defined, we form the time derivative 

V = I*® + Pi + 2^^ i 

which together with (3) can be written as 

V = e"^ (A^ P +. P A ) e + 2(6"^ P f + 6). (7) 

““III — — -i- -i_ 

Using the well-known Liapunov theorem [7], we obtain for any positive-definite 

symmetric Q, and any stability matrix A^, a positive-definite symmetric P 

as a unique solution to the equation 

-Q = P + P A . (8) 

n m 

T T 

The solution proposed by Parks [1] was to select ^(t) so that e P f. + ^ 
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Thereupon it follows fron Liapunov's direct method that V is a Liapunov 

function, and that the equilibrium at ^ ^ is asymptotically stable. This 

result does not, however, allow for certain imperfections and disturbance. 

Hence in order to define the problem which concerns us here, we write 

^2 * require that ^ be chosen so that 

P fj^ + = 0. (9) 

Then (7) together with (8) and (9) becomes 

V = -e'^Qe + 26*^ P (10) 

Here is that part of f^ which by choice or of necessity is not nullified by 
It should be noted that V is now Indefinite because the sign of 
known. To find a bound on | |^| | , a spherical region in e^ space must be 
determined such that V < 0 for ( |^1 1 ^ The derivation for vjill be 

deferred to the following section. 

At this point it is advantageous to obtain an explicit form for £2 and 

For this purpose (10) will be rewritten taking into account that f^ contains 

only one non-zero element, f . 

n 

Thus with f = f - + f „ , (10) becomes 
n nl nz 

V = -e^qe + 2Yf^2> 

n 

in which Y = I nnd, according to (9), ^ has been chosen so that 

i=l ^ 


Yf^l + = 0. 

From (4) is is seen that 
n 

f = y 6 .X. +6r-bu-d . 
n ni 1 n n n 

If the control input is now written as 

n 

= y k.x. + k . -r 
i X n+1 


( 12 ) 


(13) 


u 


(14) 
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then (13) together with (14) can be written as 
n 

f = y (6 . - b k.) X. + (<S - b k , - ) r - d . 

n ' ni nix n n n+1 n 


(15) 


We shall now assume that certain elements of the set of gains [k^, 
are adaptive, and that the remainder are identically zero. In order to satisfy 
(12) , all the terms of f^ in (15) which contain adaptive gains are used to 
comprise follows that (12) is valid if for each k^ ^ 0 there is 


an element of ^ such that 




(16) 


This result will be obtained if 




<j), = -A.b k, 
X ini 


(17) 


where is an arbitrary non-zero constant, and 


‘‘l ■ *l''/*l'’n' 

For the case in which is adaptive it follows similarly that 


(18) 


k . = ry/X ,.b . 

n+1 n+1 n 


(19) 


(14), (18), (19) represent the adaptive control law as derived in [1]. The 


point of departure from previous work is found in (10) wherein an additional 
tern due to appears in the expression for V. 


IV. Derivation of an Error Boxmd 

As a consequence of imperfect adaption and disturbance, boundedness rather 
than asymptotic stability needs to be investigated. Towards this end, a 
qp herical region is to be found such that V < 0 for ||e|| > We are now 

able to write an explicit form for (11) . Recognizing that f ^2 contains those 


1 
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paraneters which were not Identified with adaptive gains, (11) becomes 


V = -e Qe + 2Y 


ta 


■ D 

y 6x. +6r-d 
_ ^ ni i n n _ 


( 20 ) 


Here ^ signifies the sura of n terras, but not necessarily in a sequence of 
successive integers, Thus if ni=0, there is an adaptive for every state x^. 
If in addition there is an adaptive gain ^ , then the terra 6^r does not 

appear in f^^* is noted that d^ must always appear in f^ 2 * 

• T 

It is clear that V < 0 if ^ Qe > 2Yf^2‘ Thus, choosing Q = I and using 
= e^ - y^, a spherical region is sought such that, for 1 je^l | > 


T ^ o 
e e > 2 


n r n 

I Pj 1 ^ A - y.) + <Sr-d 

*^in i L ^ ni i ■'x n n _ 


( 21 ) 


(21) will be satisfied if it is required that 




r D 

n 



T r 

e e > 2 y p. 
— — ^ xn 

i=l 

l^il 

[I 

+ 1 l«„il 

'< 

H* 

+ 

6 1 jr| + 
' n' ' ' 

Denoting max |g^| = 

l%l : 

, and max 
i.t 

|h'j^(t)|* Ihj^l a 

stronger 

condition 

(22) is given by 







, . T 

f m 

m 1 


- ^ 

e e > 2p 1 i 

i=l 

le^l [i 

l“nll ( 1 

|e^l + 1 ly^l 

4 

+ 1 6 1 1 
' n' ' 

!rl + Id^l 


n 

( 22 ) 


. (23) 


Using the inequalities [8] 


ra 


n 


I le^l 1 (n I e ) 
i=l i=l ^ 


2 a/2 


? I I . ^ V 2 a/2 . , V 2a/2 

Z le.l 1 Z e ) 1 (n z 

i=l ^ 

a stronger inequality than (23) can be written as 
2 


Id I > 2^ I [ej I |6^^1 (v^ I |ej I + m |y^|) + 1 6j |i| + |d^ 


. (24) 


If 


(1 - 2^/^ ^ ° 


(25) 
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it follows that V < 0 for | |ej | > where 


R = 
e 


2v^p (q 6,1 y + 6 r + d ) 

*^in ni' -^1 n n 


1 - 2 /tun |6^^| 


( 26 ) 


The condition (25) is a test for the existence of a region < “, and is a 

sufficiency condition for stability (boundedness) of the tracking error if, 

T 

as has been assumed, ^ P ^ in (5) is positive definite. Using the notion that 

T 

an ultimate bound for ^ must be determined by a contour ^ P ^ = constant 
circumscribing the sphere of radius R^, where P is determined by (8) with 
Q = I, it can le shcvm that e must ultimately reside within a sphere of radius 

1/2 

R (27) 

e 

where max, nln eigenvalues, respectively, of the P matrix. 

Since P is positive definite, X /X , is a finite positive real number, and 

max min 

R' is finite if R is finite, 
e e 

If there is complete adaption (m=0, and 6^r does not appear in ^n2^ and if 
there is no disturbance (d^ = 0), then it is seen that R^ = 0, and the system 
is asymptotically stable in e space. 

Example of Incomplete Adaption 

A case example is introduced here to illustrate how the bound R^ can be 
applied to a meaningful problem. Suppose that the model is defined by the 
transfer function 


R' = 
e 


max 


R 


(s) = 


(s + 1)2 (|q + 1) 


10 

s^ + I2s^ + 21s + 


10 


(28) 
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and the plant is defined according to 



( + 1 )^ (^ + 1 ) 
(^i 0^2 


Cl), 

12 


3 2 2 2 

s + ( 2 oJ^ + ( 1 ) 2)3 + (cDj^ + 2 (i)^(i) 2 )s + 


( 29 ) 


The parameters ( 1)2 are assumed to lie within the ranges of 0.5 < o)^ < 1.5, 

9 < CO 2 < 11. It is clear that the coefficient (2a)j^ + has a small 
percentage variation, although can vary by ± 50%. Hence it is reasonable 
to require that and k 2 should be adaptive gains, but to assume that k^ = 0 
will not lead to excessive tracking errors. In this example the plant 
parameters were chosen to have the extreme values co^^ = 1.5, ( 1)2 = 11. 

Assuming phase variable form, it follows In (2) that 


A = 
m 


0 10 
0 0 1 
-10 -21 -12 


W'th Q = I in ( 8 ), the solution for P yields coefficients p^^^ “ 0.05, p 


23 


0.077, p ^2 = 0.048, to be used in generating y ^i 3 ^i (1®)* 

in calculating the bound in (27), the eigenvalues of P are determined to be 
\ = 0,045, 0.926, 2.507. The adaptive gains were chosen to be 


kj^ = lOyx^, k2 = lOy'x.^^ k^ = 0 

An input signal r(t) was chosen to be a square wave of unit arplitude and 1/2 
sec. period. The disturbance in this case v/as zero (^ = jO) . 

The simulation results shown in Figure 2 portray the resulting errors 
in e^(t), 62 ( 1 ), e^(t). The untimate bound on | l^l | for this case is seen 
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to be R' (actual) = l.g. 
e 

The computed bound is found from (26) , 
(27) reduces to 


R = 
e 


fa lyj) 

1 - 


(27). 


Thus in this example 


with n = 3, n = 1, p, = p„^ = 0.077. Since 6 , = 

*^in *^23 ni 

deviation term included in f » , the value for 1 6 . 

n 2 ’ ' ni 


6^2 only parameter- 

in this case becomes 


l\ll = l«33l ' Il2-2»i-“2l 

Also |y^l = max y 2 (t) =2. It follows that = 2.29 and R^ = 17 . 

As is to be expected the computed bound offers a conservative estimate. 


V. Effect of Disturbance upon Stability 

Although the results of the previous section guarantee a bound on the 
tracking error, the indefiniteness of V in (10) due to the presence of ^2 
means that stability in terns of ^ is no longer assured, even though e is 

bounded. To illustrate this point, consider the case in which e and ^ are 

scalars. In Figure 3a is depicted the solution which reaults if ^2 “ 9.' Here 
V is semidefinite, and (f) is bounded. In Figure 3b is depicted the case in 
which ^2 is nonzero. Although e is bounded by R^, it is possible that 

as shown. From (17) it is clear that |k^(->«> if , a condition which is 

unacceptable. 

We shall consider the case of complete adaption but in which a disturbance 

is present. The equation (1) for the plant can then be written as 
n 


X =y K.x,+cr+d 
n i X n n 


where K. = a , + k. , and from (18) 
X ni i 


(30) 


n 

c I 

j=l 


P. e.x. 
jn 3 i 


( 31 ) 
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where c Is an arbitrary positive constant. 

Since the purpose here is simply to demonstrate that Instability can 
occur, we shall arbitrarily choose r(t) and d(t) as step functions. With 
these inputs we see that y-*- const as t-^, and that stability of (30), (31) 
will result if e^ and K->0 (i.e. k->0) . But K-+0 requires either or 
The solution is unacceptable since this requires ( |k| |-x». Therefore we 
examine the conditions which must prevail so that e;^0. The equilibrium 
condition for (30) states that 

X = 0 = K,x, + c r + d . (32) 

n 1 1 n n 

However a necessary condition for stability at the equilibrium is that all 
are negative. Therefore it is required from (32) that 

(c r + d ) 

lim K — < 0 (33) 

1 ’‘i 

Let us consider the case in which the model output converges to r, i.e. 
y-*-r as t-x». Then with c^ ■" 1, (33) yields 

d 

— > -1. (34) 

r 

If the sign of d^ is not known, (34) can be satisfied if jr] > 

|r| < |d^| instability can result, in the sense that] 1^(34) is not 

satisfied. This result is significant in that it illustrates the danger of 

oversimplification when analyzing systems of a complex nature. 

Since the results stated above were based on the assumption of constant 

values for r and d^, it is worthwhile investigating the behavior of the 

system with disturbance when r(t) is a time varying function. Then the 

condition k = 0 will no longer represent an asymptotic solution to (31) , 

and it is possible that stability will no longer depend upon the amplitudes 

of d and r. 
n 
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Exapple the Adaption with Disturbance 

The system used in the previous example, as defined by (28), (29) will 
be used in this exatiple to show that the stability problem mentioned above 
can be avoided if the input sional is time varying. In this case there 
was complete adaption, i.e. all were adaptive, d^ was a step function, 

and the signal r(t) was a square wave. The results in Figure 4 show the 
variations in k^^(t) for various amplitudes of disturbance, and various 
square-wave frequencies. The observation is made that kj^(t) is always 
bounded and that the deviations in k^(t) become progressively smaller in 
proportion to the frequency of the square wave input. 

The stability problem discussed in this section is considered to be 
important mainly because it has been ignored. Results [9] have been reported 
recently which circumvent this problem by constructing V so that it is negative 
definite in e and For this case it is clear that both ^ and ^ will be 
bounded in presence of disturbance, and for this reason [9] is an important 
contribution. 

VI . Conclusions 

The main purpose of this investigation has been to bring to attention 
the fact that the synthesis of adaptive-control systems has often been 
discussed in the framework of idealizations which cay represent over simplifica- 
tions. A condition for boundedness of the tracking error has been derived 
for the case in which incomplete adaption and disturbance are present. However 
when using Parks 'design it is shown that instability of the adaptive 
gains can result due to the presence of disturbance. The theory has been 
applied to a non-trivial example in order to illustrate the concepts involved. 
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Tracking Errors with Incomplete Adaption 










(a) Example of Stability 



(b) Example of Instability 
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Figure 4 

Adaptive-Gain Variation with Disturbance 
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